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The electronic structure, magnetic moment, and volume collapse of MnO under pressure are 
obtained from four different correlated band theory methods; local density approximation + Hub- 
bard U (LDA+U), pseudopotential self-interaction correction (pseudo-SIC), the hybrid functional 
(combined local exchange plus Hartree-Fock exchange), and the local spin density SIC (SIC-LSD) 
method. Each method treats correlation among the five Mn 3d orbitals (per spin), including their 
hybridization with three O 2p orbitals in the valence bands and their changes with pressure. The 
focus is on comparison of the methods for rocksalt MnO (neglecting the observed transition to the 
NiAs structure in the 90-100 GPa range). Each method predicts a first-order volume collapse, but 
with variation in the predicted volume and critical pressure. Accompanying the volume collapse is a 
moment collapse, which for all methods is from high-spin to low-spin (| — » i), not to nonmagnetic 
as the simplest scenario would have. The specific manner in which the transition occurs varies 
considerably among the methods: pseudo-SIC and SIC-LSD give insulator-to-metal, while LDA+U 
gives insulator-to-insulator and the hybrid method gives an insulator-to-semimetal transition. Pro- 
jected densities of states above and below the transition are presented for each of the methods and 
used to analyze the character of each transition. In some cases the rhombohedral symmetry of the 
antiferromagnetically ordered phase clearly influences the character of the transition. 

PACS numbers: 64.30.+t,75.10.Lp,71.10.-w,71.20.-b 



I. INTRODUCTION 

For fifty years the metal-insulator transition has 
been one of the central themespj of condensed mat- 
ter physics. The type we address here does not in- 
volve spatial disorder nor change of the number of 
charge carriers per cell; the competing tendencies 
arise solely from the kinetic and potential energies in 
the Hamiltonian, favoring itineracy and localization 
respectively, and the many real-material complexi- 
ties that arise. The classic categorization is that of 
the Mott transition, treated in its most basic form 
with the single-band Hubbard model. Much has 
been learned about this model, but there are very 
few physical systems that are modeled faithfully by 
such a model. Real materials involve multiorbital 
atoms and thus extra internal degrees of freedom, 
and an environment that is often very active and 
may even react to the configuration of active sites. 

MnO is a transition metal monoxide (TMO) with 
open 3d shell that qualifies as one of the simpler 
realizations of a prototypical, but real, Mott insu- 
lator. It is, certainly, a multiorbital system with 



the accompanying complexities, but the half-filled 
3d bands lead to a spherical, spin-only moment at 
ambient pressure. Applying pressure to such a sys- 
tem leads to a number of possibilities, including 
insulator-metal transition, moment reduction, vol- 
ume collapse if a first-order transition (electronic 
phase change) occurs, and any of these may be ac- 
companied by a structural phase transition, that is, 
a change in crystal symmetry. The 3d band width 
W of such a Mott insulator is very susceptible to ap- 
plied pressure, and is one of the main determining 
factors of the strength of correlation effects. 

The "closed subshell" aspect makes MnO an atyp- 
ical 3d monoxide, as shown for example by Saito 
et at, who compiled|2j effective parameters for this 
system from spectroscopic information. An effective 
intra-atomic Coulomb repulsion energy as defined 
by them, for example, is roughly twice as large as 
for the other 3d monoxides. The complexity of this 
compound can be considered in terms of the energy 
scales that are involved in the electronic structure 
and magnetism of these oxides. These include the 
3d bandwidth W, an intra-atomic Coulomb repul- 



sion strength U, an intra-atomic d — d exchange en- 
ergy (Hund's rule J, or exchange splitting A ex ), the 
crystal field splitting A c f = e e —£t 2 , an d the charge 
transfer energy A ct = Sd — £ P (the difference in mean 
Mn 3d and O 2p site energies). In the magnetically 
ordered antiferromagnetic (AFM) state, there are 
further symmetry lowering and ligand field subsplit- 
tings involving Sd — 2p hybridization. All of these 
scales change as the volume changes, making the 
pressure-driven Mott transition a challenging phe- 
nomenon to describe. 
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FIG. 1: (Color online) A schematic (conceptual) P-T 
phase diagram for MnO based on recent high pressure 
work at Lawrence Livermore National Laboratory. 0, 0] 
The thick phase line signifies the first-order Mott tran- 
sition which simultaneously accompanies the loss of Mn 
magnetic moment, a large volume collapse, and metal- 
ization. This transition should end at a critical point 
(solid circle). The gray fan above the critical point sig- 
nifies a region of crossover to metallic behavior at high 
temperature. The star denotes earlier shock data,|| and 
the diamond marks the ambient pressure Neel tempera- 
ture. Only the distorted Bl (dBl) phase is magnetically 
ordered. 

Although the objective of the current paper is to 
compare methods within the fee (rocksalt) phase, it 
is useful first to recount what is known about the 
Mott transition at this time. The current experi- 
mental information, mostly at room temperature, on 
the behavior of MnO under pressure is summarized 
in Fig. ^ Resistance measurements |3| provided the 
first evidence of the Mott transition in MnO near 100 
GPa. Recent x-ray diffraction and emission spec- 
troscopy measurements of the crystal structure and 
magnetic moment by Yoo et aL[4| have clarified the 
behavior. Around 90 GPa there is a structural trans- 
formation from the distorted Bl (rocksalt) phase to 
the B8 (NiAs) structure. This structure change is 
followed at 105 GPa by the Mott transition, consist- 
ing of a simultaneous volume collapse and moment 
collapse signifying a qualitative change in the elec- 
tronic structure of the compound. 



On the theoretical side, little is known about how 
the Mott transition occurs in a real multiband TMO 
in spite of the extensive studies of the Mott transi- 
tion in the single-band Hubbard model, which has 
a simple spin-half moment at strong coupling and 
half-filling. The numerous energy scales listed above, 
and the S=| moment on Mn arising from the five 
3d electrons, allow many possibilities for how the 
moment might disintegrate as the effective repulsion 
decreases. The high pressure limit is clear: a non- 
magnetic 3d—2p band metal in which kinetic energy 
overwhelms potential energy. This is the competi- 
tion studied in the (simplified) Hubbard model. The 
multiband nature has attracted little attention until 
recently, when for example the question of possible 
orbital-selective Mott transitions |y, L3 have aroused 
interest. One can imagine one scenario of a cas- 
cade of moment reductions S = § — > § — ► ^ before 
complete destruction of magnetism, as electrons use 
their freedom to flip spins (as some competing en- 
ergy overcomes Hund's rule, for example). In such a 
scenario there is the question of which orbital flips 
its spin at each spinflip, which involves a question of 
orbital selection and ordering. At each flip the sys- 
tem loses exchange (potential) energy while gaining 
kinetic energy (or correlation energy through 'sin- 
glet' formation). The manner in which kinetic en- 
ergy changes is difficult to estimate because subband 
involvement means that there is no longer a single 
bandwidth W that is involved. The increasing hy- 
bridization with O 2p states under pressure strongly 
affects the kinetic energy, directly and through su- 
perexchange (a kinetic energy effect). 

It is known that conventional band theory (local 
density approximation [LDA]) that does so well for 
so many materials gives poor results for 3d monox- 
ides in many respects, and some predictions are qual- 
itatively incorrect (viz. no band gap when there 
should be a large gap of several eV) . Thus even at the 
density functional level (ground state energy, den- 
sity, and magnetization) some correlated approach 
is required. In the past fifteen years several ap- 
proaches, which we refer to as correlated band theo- 
ries, have been put forward, and each has had its suc- 
cesses in providing an improved description of some 
aspects of correlated TMOs. Although commonly 
called mean-field approaches with which they share 
many similarities, they are not mean-field treat- 
ments of any many-body Hamiltonian. Rather, they 
are energy functionals based on the complete many- 
body Hamiltonian, which must then be approxi- 
mated due to limited knowledge of the exchange- 
correlation functional. 

In this paper, we provide a close comparison 
of certain results from four such methods: full 
potential LDA+U, the hybrid exchange functional 
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(HSE) approach, the self-interaction-corrected lo- 
cal spin-density method (SIC-LSD), and a nonlo- 
cal pseudopotential-like variation of SIC (pseudo- 
SIC). Our main focus is to compare the predicted 
changes in energy, moment, and volume within the 
AFMII rocksalt phase of MnO. To keep the com- 
parison manageable we confine our attention to the 
rocksalt phase, since our emphasis is on compari- 
son of methods and not yet the ultimate but more 
daunting task of modeling structural changes that 
may precede, or accompany, the Mott transition. 



II. STRUCTURE AND SYMMETRY 

Rocksalt MnO has an experimental equilibrium 
lattice constant a D =4.45 A (cubic cell volume 
V =88.1 A 3 ). Density functional theory, like 
Hartree-Fock theory, deals in its most straightfor- 
ward form with ground state properties, i.e. zero 
temperature. The ground state is known to be the 
AFMII phase in which < 111 > layers have spins 
aligned, and successive layers are antiparallel. The 
resulting symmetry is rhombohedral, with Mn| and 
MnJ. being distinct sites (although related through a 
translation + spin- flip operation). Thus, while most 
of the lore about transition metal monoxides is based 
upon cubic symmetry of the Mn (and O) ion, in the 
ordered state the electronic symmetry is reduced. It 
is obvious that individual wavefunctions will be im- 
pacted by this symmetry, viz. fourfold symmetry 
around the cubic axes is lost. It has been empha- 
sized by Massidda et al. 8] that zone-integrated, and 
even spin-integrated, quantities show the effects of 
this symmetry lowering; for example, Born effective 
charges lose their cubic symmetry. Since this issue 
arises in the interpretation of our results, we provide 
some background here. 

In cubic symmetry the Mn 3d states split into the 
irreducible representations (irreps) denoted by t2 g 
and e g . Rhombohedral site symmetry results in the 
three irreducible representations a g , e g> i, and e 9i 2, 
the latter two being two fold degenerate. The coor- 
dinate rotation from cubic to rhombohedral (super- 
script c and r respectively) is, with a specific choice 
for the orientation of the x and y axes in the rhom- 
bohedral system, 
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Applying this rotation of coordinates gives the 3d 
orbitals in the rhombohedral frame in terms of those 
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In rhombohedral coordinates it is useful to cate- 
gorize the 3d orbitals in terms of their orbital an- 
gular momentum projections along the rhombohe- 

±1; 

mi — ±2. It is easy to see that 
mi\ specifies groups of states that only transform 
into combinations of themselves under trigonal point 
group operations. 

Note that the unique a g symmetry state in rhom- 
bohedral coordinates is the fully symmetric combi- 
nation of the cubic t2 g states. The other two irreps 
are both e g doublets. While \mg\ — 1 and \mg\ = 2 
form representations of these irreps, if there are com- 
ponents of the crystal field that are not diagonal in 
the L=(2,m^) basis, these states will mix. Then each 
of the resulting (orthonormal) irreps e ffl i, and e 9i 2 
will contain both \me\ — 1 and \mi\ — 2 components. 
Such mixing does occur in MnO and complicates the 
symmetry characterization of the 3d states. 



III. METHODS 

A. LDA Calculations 

For LDA band structure plot (Fig. 0) we used 
version 5.20 of the full-potential local orbital band 
structure method (FPLOJ3)- Relativistic effects 
were incorporated on a scalar-relativistic level. We 
used a single numerical basis set for the core states 
(Mn Is2s2p and O Is) and a double numerical ba- 
sis set for the valence sector including two 4s and 
3d radial functions, and one 4p radial function, 
for Mn, and two 2s and 2p radial functions, and 
one 3d radial function, for O. The semi core states 
(Mn 3s3p) are treated as valence states with a sin- 
gle numerical radial function per nl-sheW. The lo- 
cal density exchan ge-c orrelation functional PW92 of 
Perdew and Wang|10| was used. 
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B. LDA+U Method 

The LDA+U approach of including correlation ef- 
fects is to (1) identify the correlated orbital, 3d in 
this case, (2) augment the LDA energy functional 
with a Hubbard-like term (Coulomb repulsion U) 
and Hund's (exchange J) energy between like spins, 

(3) subtract off a spin-dependent average of this in- 
teraction energy to keep from double-counting repul- 
sions (once in LDA fashion, once in this U term), and 

(4) include the correlated orbital occupation num- 
bers in the self-consistency procedure, which leads 
to an orbital-dependent Hartree-Fock-like potential 
acting on the correlated orbitals. The addition to 
the energy functional has the schematic form 

E v = - ^2 (U - JS a a')[n mC rn m 'a' - n CT rv].(6) 

rcia^m' a' 

We actually use the coordinate-system indepen- 
dent form of LSDA+U [U El El implemented 
in FPLO,|9j which leads to four m indices on U 
and J which for simplicity have not been displayed 
(nor has the full off-diagonal form of the occupation 
matrices n mm i a ). This treatment of the on-site 
interactions U and J incorporates on-site corre- 
lation effects in the Mn 3d-shell. We have used 
the so called 'atomic-limit' (strong local moment) 
form of the double-counting correction, the last 
term in Eq. [BJ This form is appropriate for the 
high-spin state, but it is less obviously so for the 
low-spin state that is found at reduced volumes. 
The Slater parameters were chosen according to 
U = F = 5.5 eV, J = i(F 2 + F A ) = 1 eV and 
F 2 /F 4 = 8/5. 

The shape of the basis orbitals has been opti- 
mized yielding a sufficient accuracy of the total en- 
ergy over the range of geometries considered in this 
work. The k-integrals are performed via the tetrahe- 
dron method with an irreducible mesh corresponding 
to 1728 (12 3 ) points in the full Brillouin zone. 

C. SIC-LSD Method 

The SIC-LSD method addresses the unphysical 
self-interaction in the LDA treatment of localized 
states. Itinerant states, being spread over space 
without finite density in any given region, do not ex- 
perience this self-interaction within the LDA treat- 
ment. Should there be localized states, confined 
to some region and giving a finite density, they 
will suffer an unphysical self-interaction in the LSD 
method. This issue then clearly arises in the 
itinerant-localized transition in MnO and other cor- 



related systems. The basic premise of the SIC- 
LSD method is that localized electrons should ex- 
perience a different potential from that of itinerant 
electrons . [Til IT5L IT^| analogous to that of an atomic 
state whose self-interaction must be removed. Then 
electrons on the surrounding atoms are allowed to 
accommodate self-consistently. This distinction of 
localized versus itinerant state is addressed in SIC- 
LSD by extending the energy functional in the form 

occ. 

E SIC-LSD = E LSD_J2 5 SIC : (?) 

a 

€° = U[n a ]+E^ D [n a }. 

Here U[n] represent the Hartree (classical Coulomb) 
energy of a density n(r). The self-Coulomb en- 
ergy U[n a ] and self-exchange-correlation E^ D [n a ] 
energies are subtracted off for each localized state 
ip a with density n a . Whether states are local- 
ized or not (with non-zero, respectively zero self- 
interaction) is determined by minimization of this 
functional, allowing localized as well as itinerant 
states ip a . Since the correction vanishes for itiner- 
ant states, the sum finally includes only the self- 
consistently localized states. The localized and itin- 
erant states are expanded in the same basis set, 
and minimization becomes a process of optimizing 
the coefficients in the expansion of the states (as 
other band structure methods do, except that Bloch 
character is imposed in other methods). The im- 
plementation of Temmerman and collaborators^^ 
used here incorporates the atomic-sphere approxi- 
mation (ASA) of the linear muffin-tin orbital elec- 
tronic structure method 18] (LMTO) in the tight- 
binding representation. [19| Further details can be 
found in Ref. |U 

D. pseudo-SIC Method 

The large computing requirements (compared to 
LDA) of the SIC-LSD method, even for materials 
with small unit cells, has led to an alternative ap- 
proach [2lJ , in which the self- interaction part of the 
Kohn-Sham potential is approximated by a non- 
local, atomic-like contribution included within the 
pseudopotential construction. The original imple- 
mentation of this scheme has given important im- 
provements over LSDA results for non-magnetic II- 
VI and III-V semiconductors, but was not applica- 
ble to metals or to magnetic and highly-correlated 
systems where there is a coexistence of strongly lo- 
calized and hybridized electron charges. 

The pseudopotential self-interaction corrected 
calculations presented here were performed using 
the recently developed "pseudo-SIC" method of 
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Filippetti and Spaldin [22]. This pseudo-SIC 
approach represents a compromise between the 
fully self-consistent implementations of Svane et 
al. 0, IsS and the alternative method of Vogel 
et al. |2l|, in that the SIC calculated for the 
atom (as in Ref. 0|) is scaled by the electron 
occupation numbers calculated sclf-consistently 
within the crystal environment. This allows the SIC 
coming from localized, hybridized, or completely 
itinerant electrons to be discriminated, and permits 
the treatment of metallic as well as insulating 
compounds, with minimal computational overhead 
beyond the LSDA. In this pseudo-SIC procedure, 
the orbital SIC potential is taken from the isolated 
neutral atom and included in the crystal potential 
in terms of a nonlocal projector, similar in form 
to the nonlocal part of the pseudopotential. The 
Bloch wave functions are projected onto the ba- 
sis of the pseudo-atomic orbitals, then, for each 
projection, the potential acting on the Bloch state 
is corrected by an amount corresponding to the 
atomic SIC potential. Note that, within this for- 
malism, a physically meaningful energy functional 
which is related to the Kohn-Sham equations by a 
variational principle is not available. However, a 
suitable expression for the total energy functional 
was formulated in Ref. and shown to yield 

structures in good agreement with experiment. We 
use this functional here to calculate the bondlength 
dependence of the total energy. We have used 
ultra-soft pseudopotentials with an energy cutoff of 
35 Ry. An 8 3 Monkhorst-Pack grid was used for 
fc-point sampling. The low-spin and high-spin solu- 
tions were obtained by setting initial magnetization 
to 5/ib or lfiB respectively. 



E. Hybrid Functional Method 

The hybrid-exchange DFT approximation mixes a 
fraction of the exact, non-local, exchange interaction 
(which uses the Hartree-Fock [HF] expression) with 
the local, or semi-local, exchange energy of the LDA 
or the generalized gradient approximation (GGA). 
The PBEO functional takes the form: 



E xc = aE* F + (1 - a)E™ m + E, 



PBE 



pPBE 



(8) 



where E BF and E PBE are the exchange and the 
PBE GGA functionals. The mixing parameter a = 
1/4 was determined via perturbation theory [21] and 
E PBE is the PBE correlation energy. 

In this work we use the hybrid method recently de- 
veloped by Heyd, Scuseria and Enzerhof (HSE)).[24| 
It is based upon the PBEO functional, but employs 



a screened, short-range (SR) Hartree-Fock (HF) ex- 
change instead of the full exact exchange, which re- 
sults in a more efficient evaluation for small band 
gap systems. In this approach, the Coulomb oper- 
ator is split into short-range (SR) and long-range 
(LR) components respectively 



1 1 — erf (lot) erf(u>r) 
r r r 



(9) 



where w is a parameter that can be adjusted for 
numerical or formal convenience. 

The expression for the HSE exchange-correlation 
energy is 



E Tr — clE 



HF,SR 



-E 



PBE,LR 



(w) + (1 - a)E 

H + e, 



PBE.SR 



M(io) 



PBE 



where F, bf,sr (uj) is the SR HF exchange com- 
puted for the SR part of the Coulomb potential, 
E PBE ' SR (uj) and E pbe ' lr (uj) are the SR and the 
LR components of the PBE exchange, respectively. 
The cited papers should be consulted for further de- 
tails. The HSE functional has been found to yield 
results in good agreement with exp eriment for a wide 
range of solids and molecules |25l l2fi| . 

This functional is implemented in the develop- 
ment version of the Gaussian quantum chemi stry 
package [H. We use the Towler basis [H H3, 
of Gaussian functions for our basis set. It consists 
of a [20sl2p5d/5s4p2d] contraction for Mn and a 
[14s6p/4s3p] basis for O, optimized for HF stud- 
ies on MnO. A pruned grid for numerical integra- 
tion with 99 radial shells and 590 angular points per 
shell was used. The fc-space was sampled with a 16 3 
mesh. Low spin and high spin antifcrromagnetic ini- 
tial guesses were obtained using the crystal field ap- 
proach by patching the density matrix obtained from 
diagonalization of the Harris functional j^] with the 
density matrices obtained from calculation on ions 
in the appropriate ligand field |33j, |3J] . 



IV. PREVIOUS ELECTRONIC 
STRUCTURE STUDIES 

The origin, and the proper description, of the 
moments and the band gaps in transition metal 
monoxides have a long history. The earliest ques- 
tion centered on the connection between the antifcr- 
romagnetic (AFM) order and the insulating behav- 
ior. Slater's band picture [35| could account in a one- 
electron manner for a gap arising from AFM order, 
whereas Mott's picture of correlation-induced insu- 
lating behavior 36] was a many-body viewpoint with 
insulating behavior not connected to the magnetic 
order. The proper general picture in these monox- 
ides arose from studies of transport above the Neel 
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temperature and with introduction of defects, giving 
them the designation as Mott insulators. 

Much progress on the understanding of MnO 
and the other monoxides came from early studies 
using LDA. While understanding that LDA does 
not address the strong correlation aspect of the 
electronic structure, Mattheiss[37j and Terakura et 
al. 38] quantified the degree and effects of 3d — 2p 
interactions, and pointed out the strong effect of 
magnetic ordering on the band structure. More re- 
cently, Pask and collaborators 39] have studied the 
structural properties, and the rhombohedral distor- 
tion, with LDA and GGA approximations. The sym- 
metry lowering and resulting structure is described 
well, and in addition they found that AFM order- 
ing results in significant charge anisotropy. Effects 
of AFM order were further probed by Posternak et 
al. by calculating and analyzing maximally localized 
Wannier functions for the occupied states. |40| 

The application of correlation corrections in MnO 
already has a colorful history. The first work, by 
Svane and Gunnarsson|23j and by Szotek et aL,|4lj 
was in the application of the SIC-LSD method. The 
former pair correctly obtained that MnO, FeO, CoO, 
NiO, and CuO are AFM insulators, while VO is a 
metal. They calculated a gap of 4 eV for MnO. 
Szotek et al. used a fairly different implementa- 
tion of the SIC-LSD approach but find a similar 
gap (3.6 eV). Their 3d states lay about 6 eV below 
the center of the 2p bands, although hybridization 
was still clearly present. In this same time frame, 
Anisimov, Zaanen, and Andersen introduced|42|| the 
LDA+U method with application to the transition 
metal monoxides. They obtained a band gap of 3.5 
eV but few other results on MnO were reported. 

Kotani implemented[H |H the "DFT exact- 
exchange" method of Talman and Shadwick^fj to 
crystal calculations. This method consists of tak- 
ing the Fock expression for the exchange energy in 
the DFT functional, then performing a Kohn-Sham 
solution (minimization), giving a local exchange po- 
tential ("optimized effective potential"). In Kotani's 
results for MnO, the Mn e g and t2 g bands form very 
narrow (almost atomic-like) bands between the oc- 
cupied O 2p bands and the conduction bands. Taka- 
hashi and Igarashi|47j proposed starting from the 
Hartree-Fock exchange and adding correlation from 
a local, three-body scattering viewpoint. Their cor- 
rections were built on a parametrized tight-binding 
representation, and they obtained small self-energy 
corrections for MnO, much smaller than they ob- 
tained for the other transition metal monoxides. 

The effective potential approach used by Kotani 
was extended by Solovyev and Terakura 48] in an 
unconventional way. They obtained an effective po- 
tential using the criterion that it had to reproduce 



the spin- wave spectrum, i.e. that it had to describe 
the magnetic interactions correctly. They found 
clear differences when comparing to the LDA+U 
and the optimized effective potential results, and dis- 
cussed limitations of the one-electron band method 
itself. 

More recently, Savrasov and Kotliar applied a dy- 
namical extension^! of the LDA+U method (dy- 
namical mean field theory) to MnO and NiO. Be- 
ing a self-energy method, this is not really a cor- 
related band theory. For the properties they cal- 
culated (band gap, effective charges, dielectric con- 
stant, optic phonon frequencies) the dynamical re- 
sults are similar to the LDA+U results and differ 
considerably from LDA values. 

Even though hybrid-exchange DFT applications 
to solids are still in their infancy, there have been two 
previous studies of MnO. The first, by Bredow and 
Gerson,[50j utilized the B3LYP hybrid functional. 
Unlike the LDA and GGA, they found B3LYP pro- 
vided an excellent band gap for MnO. More recently, 
Franchini et al. have examined MnO in more de- 
tail using the PBEO approximation. |51j They also 
found a gap, lattice constant and density of states 
in quite good agreement with experiment. In par- 
ticular, the distorted dBl rhombohedral structure 
was determined to be the minimum energy geometry, 
in agreement with experiment. Neither the B3LYP 
nor the PBEO approximation can be applied to the 
metallic side of the transition of interest here. For 
that, we must turn to the screened hybrid-exchange 
of HSE. 

Therefore, while there has been thorough LDA 
studies of MnO and a variety of approaches to treat- 
ment of the correlation problem, nearly all of these 
have considered only ambient pressure or small vari- 
ations of volume near zero pressure. The work de- 
scribed in the following sections focuses on testing 
the four different correlated band methods from am- 
bient conditions to high pressures, through the vol- 
ume collapse regime, to see whether some basic foun- 
dation can be laid for the understanding and theoret- 
ical description of pressure-driven Mott transitions 
in real materials. 



V. RESULTS 

Our principal results revolve around the first- 
order transition, high— dow volume which is also 
high— dow moment in nature. For convenience we 
use from here on the specific volume v = V/Vq, the 
volume referenced to the experimental zero-pressure 
volume. The equations of state have been fitted for 
both high volume and low volume phases for each 
computational method, and the resulting constants 



G 





FPLO 


Gaussian 


PW-USSP 


LMTO-ASA 


a/a a 


0.97 


0.97 


0.96 


0.97 


B (GPa) 


170 


196 


169 


205 


Gap (eV) 


0.72 


1.13 


0.92 


1.04 


Moment (/xs) 


4.52 


4.53 


4.42 


4.42 



TABLE I: Local density approximation results obtained 
from the four codes used in this work, for the equilib- 
rium lattice constant relative to the experimental value 
(a/a ), the bulk modulus B, the energy gap, and Mn 
moment. The experimental values for B fall in the 
range 142-160 GPa (see caption to Table Ell ■ AH cal- 
culations are carried out for the antiferromagnetic II 
phase of MnO in the cubic (NaCl) structure. Differ- 
ences are due to (1) differing exchange-correlation func- 
tions that are used, (2) approximations made in algo- 
rithms to solve the Kohn-Sham equations and evaluate 
the energy, (3) basis set quality, and (4) for the moment, 
the operational definition of the "Mn moment" differs. 
The notation follows. FPLO: full potential local orbital 
method with Perdew-Wang (1992) exchange-correlation 
functional. Gaussian: Gaussian local orbital code, Slater 
exchange (a=2/3) and Vosko-Wilk-Nusair correlation. 
PW-USSP: planewave basis using ultrasoft pseudopoten- 
tials, with the Perdew-Zunger exchange-correlation func- 
tional. LMTO-ASA: linear muffin-tin orbital code in the 
atomic sphere approximation, using the Perdew-Zunger 
exchange-correlation functional. 
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FIG. 2: (Color online) LDA band structure of AFM 
MnO along rhombohedral symmetry lines, calculated 
with the FPLO method|?|, with horizontal line ("Fermi 
level") placed at the top of the gap. The F-T lies along 
the rhombohedral axis, while F-L lies in the basal plane. 
The O 2p bands lie in the -8 eV to -3.5 eV range, with 
the majority Mn 3d bands just above (-3 eV to -1 eV). 
The five minority 3d bands are just above the gap. Note 
the small mass, free-electron-like band that lies below 
the unoccupied 3d bands at the F point. 



are presented and analyzed below. 



A. Baseline: LDA Bands and Equation of 
State 

In Table[I]we illustrate the magnitude of variation 
of three properties calculated from the four codes 
used here, as an indication of what size of differ- 
ences should be given meaning in properties pre- 
sented below. Before the method-specific beyond- 
LDA corrections are applied, there are differences in 
the bandgap (0.95±0.2 eV) due to the algorithms 
applied in the codes, to the basis quality, and be- 
cause different LDA functionals are used. It is evi- 
dent that these variations have very little effect on 
the calculated equilibrium volume and minor effect 
on the calculated Mn moment (4.47±0.05 hb)- Such 
differences will not affect the comparisons of the re- 
sults given here, at the level of precision of interest 
at this time. 

The LDA band structure of AFM MnO is shown 
in Fig. |21 as the reference point for the following 
calculations. There is a band gap of ~ 0.7 eV. The 
five bands immediately below the gap are the ma- 
jority Mn 3d bands, those lying below are the O 



2p bands. The charge transfer energy mentioned 
in the Introduction is A ct = Sd — £ P = 6 eV, and 
the exchange splitting is ~3.5 eV. It is tempting to 
interpret the 3+2 separation of occupied 3d states 
as t2 g + e g , but the rhombohedral symmetry ren- 
ders such a characterization approximate. The five 
bands above the gap are primarily the minority Mn 
3d bands. However, a free-electron-like band at T 
lies lower in energy than the 3d bands, but disperses 
upward rapidly, so over most of the zone the low- 
est conduction band is Mn 3d and the gap is 1 eV. 
The presence of the non-3d band does complicate 
the interpretation of the band gap for some of the 
correlated methods, presented below. 

The behavior of MnO under compression within 
GGA has been given earlier by Cohen, Mazin, and 
Isaak. 53J They obtained an equilibrium volume 2% 
higher, and bulk modulus 13% smaller, than mea- 
sured. Pressure studies including extensive struc- 
tural relaxation have also been provided by Fang 
et aL[54| Their structural relaxations make their 
study more relevant (within the restrictions of GGA) 
but also make comparison with our (structurally re- 
stricted) results impossible. 
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B. Energetics and Equation of State 





GGA 


LDA+U 


HSE 
exchange 


pseudo 
-SIC 


SIC 
-LSD 


Vo 


1.02 


0.93 


0.99 


1.09 


1.04 


Vh 


0.70 


0.66 


0.60 


0.86 


0.64 


Vl 


0.62 


0.61 


0.55 


0.73 


0.52 


Av 


0.08 


0.05 


0.05 


0.13 


0.12 


B h 


196 


192 


187 


138 


159 


Bi 




195 


224 


230 


67 


B' h 


3.9 


3.2 


3.3 


3.6 


3.3 


Bi 




3.6 


4.0 


3.5 


4.7 


Pc 


149 


123 


241 


56 


204 



TABLE II: Quantities obtained from fits to the Mur- 
naghan equation of state for the various functionals, ex- 
cept for the GGA column, which are taken from Ref. 
[53. vo is the calculated equilibrium volume, B and B' 
are the bulk modulus (in GPa) and its pressure deriva- 
tive, v h , vi are the calculated volumes of the high and 
low pressure phases, respectively, at the critical pressure 
P c (in GPa). An is the amount of volume collapse that 
occurs at the transition pressure P c . All volumes are 
referred to the experimental equilibrium volume. The 
experimental values are 5=142-160 GPa, B' ~ 4; see 
Zhang|52^| and references therein. 

The equation of state (EOS) energy vs. volume 
curves for the various functionals are collected in 
Fig. |3J For each correlated band method a large 
volume, high-spin state and a small volume, low-spin 
state are obtained. The analysis to obtain the first- 
order volume collapse transition was done as follows. 
For each high volume and low volume phase sep- 
arately, an EOS function Eh t iiV) was determined 
(h, i=high, low) by a fit to the Murnaghan equa- 
tion. Both fits give minima, with the most relevant 
one being for the high spin phase and being the pre- 
dicted equilibrium volume Vq . The pressure is ob- 
tained from the volume derivative of the EOS, which 
is inverted to get V(P). Equating the enthalpies 
E[V(P)] + PV{P) of the two phases gives the crit- 
ical pressure P c . The volumes at this pressure then 
give the volume collapse AV = Vh(P c ) — Vj(P c ). 

The various quantities for all four computational 
schemes are given in Table [H] along with the mi- 
correlated results of Cohen et a/.|53| Not surpris- 
ingly given the other differences that will be dis- 
cussed, there are substantial variations among the 
critical pressures and related quantities. Particu- 
larly noticeable already in the EOS is the result that 
the energy difference between the low-spin energy 
minimum and the high-spin one is ~0.2 eV for the 
LDA+U and pscudo-SIC methods, while the HSE 




V/V (a = 4.45A) V/V (a c = 4.45 A) 

FIG. 3: (Color online) The calculated total energy/MnO 
versus volumes for the various functionals for the AFMII 
rocksalt phase, referred to zero at their equilibrium vol- 
ume. The hlled symbols denote the calculated energies 
of the large volume, high spin configuration and the open 
symbols denote the calculated energies of the small vol- 
ume, low spin configuration. The continuous and dashed 
lines are the least square fitted curves to the Murnaghan 
equation of state for high and low spin configurations 
respectively. 



method gives roughly twice the energy difference 
(0.4 eV), and the SIC-LSD method gives roughly 
0.6 eV. 

We now mention other noteworthy features of the 
calculated data in Table ITTI 

(1) The predicted equilibrium volume from the 
LDA+U method is the smallest of the four meth- 
ods (u^=0.93), thus overbinding. The HSE value 
is almost indistinguishable from the observed value, 
while the SIC-LSD and pseudo-SIC methods give 
underbinding (vq 1 = 1.04, 1.09 respectively). 

(2) The pseudo-SIC method predicts the transi- 
tion to occur at a relatively small volume reduction 
(v h =0.86); the other methods give the onset of tran- 
sition at v=0. 63+0.03. 

(3) The critical pressure P c =56 GPa predicted by 
pseudo-SIC is smallest of the methods. P c in 
LDA+U (123GPa) is comparable to that of LDA; 
those of SIC-LSD and HSE are higher (204, 241 GPa, 
respectively). 

(4) The SIC-LSD method predicts a transition to a 
low volume phase that is much softer than the high 
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volume phase, a phenomenon that is extremely un- 
usual in practice but not disallowed. There are two 
possible sources of this difference: (i) in SIC-LSD 
the system becomes completely LDA-like in the low 
volume phase whereas in the other methods the 3c? 
states are still correlated, or (ii) the LMTO-ASA 
method involves approximations that pose limita- 
tions in accuracy. |55| 

(5) The values of B in the large volume phase vary 
although not anomalously so, given the differences 
discussed just above. 

(6) The values of B' in the high volume phase 
are reasonably similar across the methods: B' = 
3.4 ± 0.2. The variation in the collapsed phase is 
greater. 



C. Magnetic Moment 

The moment collapse behavior of each method is 
collected in Fig. 0] For comparison, the GGA re- 
sult presented by Cohen et aL[53j was a moment 
collapse from 3.4/is to 1.3/j.b at the volume given 
in Table [Q] On the broadest level, the predictions 
for the Mn moment show remarkable similarity. At 
low pressure, all methods of course give the high- 
spin S=| configuration of the Mn 2+ ion, with the 
local moment being reduced slightly from 5/1^ by 
3d — 2p mixing. This electronic phase persists over a 
substantial volume reduction, giving way in all cases 
to an state, not the nonmagnetic S = result 
that might naively be anticipated. Three methods 
give a stable moment very near over a range 
of volumes. The SIC-LSD method (which in the 
collapsed phase is simply LDA) is alone in giving a 
varying moment, one that reduces from 1.4 fiB at 
v = 0.68 down to 0.5 [i B at v = 0.48. 

One difference between the methods lies in how 
soon the low-spin state becomes metastable, i.e. 
when it is possible to obtain that state self- 
consistently, as opposed to when it becomes the sta- 
ble solution (which was discussed in the EOS subsec- 
tion). The state is obtained already at ambient vol- 
ume in the HSE method; the pseudo-SIC method ob- 
tains the low-spin state just below v — 0.80; for the 
LDA+U method, it was not followed above v = 0.68. 
It should be emphasized however that no concerted 
measures were taken to try to follow all solutions to 
the limit of their stability. 



D. Fundamental Band Gap 

In Fig. [3] the calculated band gap of both high- 
spin and low-spin states are shown for all methods. 
Here the behavior differs considerably between the 
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FIG. 4: (Color online) Calculated values of the mo- 
ment (shown as symbols) on each Mn site as a function 
of volume for the various functionals. All the methods 
predict a distinct collapse (first order) of magnetic mo- 
ment with decrease in volume. At large volumes, the 
high spin state with S = 5/2 (single occupancy of each 
of the 3d orbitals) is realized while the low spin state 
with S = 1/2 is favored for smaller volumes. Note: the 
computational methods calculate the 'Mn moment' in 
inequivalent ways, so differences of less than 0.1/is may 
have no significance. 



methods, in part because at certain volumes the gap 
lies between different bands for some of the meth- 
ods. At ambient pressure the pseudo-SIC and HSE 
methods obtain a gap of 3.5-4 eV, while that for 
SIC-LSD is 2.9 eV, and that of the LDA+U method 
is even lower, less than 2 eV. Experimental values 
lie in the 3.8-4.2 eV range. Referring to Fig. 2, it 
can be observed that the large volume gap depends 
on the position of the majority 3d states with re- 
spect to the free-electron-like band, i.e. it is not the 
3d — 3d Mott gap. Both of the former approaches 
show only a slight increase as pressure is applied, 
reaching a maximum around v — 0.76 where a band 
crossing results in a decreasing gap from that point. 
For pseudo-SIC, there is an almost immediate col- 
lapse to a metallic low spin state. Within HSE, MnO 
collapses to the low-spin state at a volume v = 0.55, 
coincident with the closing of the gap in the low- 
spin state. The SIC-LSD and LDA+U gaps, smaller 
initially, show a much stronger increase with pres- 
sure, and also incur the band crossover that leads to 
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FIG. 5: (Color online) Calculated band gap as a function 
of volume; the lines are simply a guide to the eye. For 
LDA+U, the band gap increases with decrease in volume 
for the high spin state, but decreases with volume in the 
low spin state. For the HSE calculations, the band gap 
monotonously decreases with volume whilst pseudo-SIC 
shows a first order jump in going from high to low spin 
state. At very low volumes, where the low spin configu- 
ration is preferred LDA+U gives a substantial gap and 
is still an insulator, while HSE and pseudo-SIC calcula- 
tions converge to a metallic solution. 



decrease of the gap (within the high-spin state). 



VI. ANALYSIS OF THE TRANSITION 

In this section we analyze the character of the 
states just above and just below the Mott transi- 
tion, as predicted by each of the methods. Due to 
the differing capabilities of the codes, the quantities 
used for analysis will not be identical in all cases. In 
the Figs. 6-9 we present for uniformity the DOS 
in the high volume phase at the equilibrium vol- 
ume (ao) and in the collapsed phase at a — 0.85ao 
(v=0.6). Note (from Table I) that this specific vol- 
ume does not correspond to any specific feature in 
the phase diagram for any method, although it lies 
in the general neighborhood of the volume at the 
collapse. Changes within the collapsed phase are 
continuous, however, so the plots at a = 0.85ao are 
representative of the collapsed phase. 
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FIG. 6: (Color online) Projected DOS onto symmetrized 
Mn 3d orbitals in the rhombohedral AFMII rocksalt 
phase using the LDA+U method. Top panel: High spin 
solution at the LDA+U equilibrium volume. Bottom 
panel: Low spin solution at 60% of the LDA+U equilib- 
rium volume. The a g orbital is the 3z 2 — r 2 oriented along 
the rhombohedral axis, other symmetries are described 
in the text. The overriding feature is the spin-reversal of 
the m = ±1 e' g orbitals between the two volumes. 



A. LDA+U 

The projected DOSs (PDOSs) in Fig. © refer to 
projections onto Mn 3d orbitals, with the z-axis be- 
ing the rhombohedral axis, the a g 3z 2 - r 2 (|m|=0) 
state; the e' g pair {xz,yz} (|m|=l); and the e g pair 
{x 2 — y 2 ,xy} (\m\=2). Because the two e g repre- 
sentations have the same symmetry, they can mix 
and the actual combinations e 9i i, e Si 2 are orthogonal 
linear combinations of e g , e' g which depend on inter- 
actions. For the LDA+U results, however, there is 
little mixing of the e g , e' g pairs. The character of the 
transition is simple to describe: the e' g pair (|m|=l 
with respect to the rhombohedral axis) simply flips 
its spin. 

This S=o state is unexpected and quite unusual. 
First, each 3d orbital is still singly occupied, veri- 
fied by plotting the charge density on the Mn ion 
and finding it just as spherical as for the high- 
spin state. Second, each 3d orbital is essentially 
fully spin-polarized, with the configuration being 
a g f e g | e g I- A plot of the spin density [5(J reveals 
the unanticipated strong anisotropy with nodal char- 
acter, characteristic of spin-up m = and |to|=2 or- 
bitals, and spin-down |m| = l orbitals (in the rhombo- 
hedral frame). Third, it makes this transition with 
essentially zero change in the gap, which is 3.5 eV. 
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FIG. 7: (Color online) Spin- and orbital-projected DOS 
from the HSE (hybrid-exchange) method. Orbitals are 
expressed in a coordinate system that is specific to the 
code, neither cubic nor rhombohedral. In the collapsed 
phase two orbitals are (nearly) doubly occupied, and the 
net moment arises from the single occupation of the or- 
bital labeled 1I4. In a cubic representation the configura- 
tion can be specified as t\ g = t% g ^ + t\ g y . At the volume 
u=0.55 of onset of the collapsed phase (see Table I), the 
gap is essentially zero. 



The band structure changes completely, however, so 
the close similarity of the gaps on either side is ac- 
cidental. 



B. HSE Method 



In the high volume phase the distribution and 
overall width of the occupied 3d states, shown in Fig. 
is similar to that of the LDA+U method (previous 
subsection). The gap is larger, as discussed earlier. 
The collapsed phase shows new characteristics. The 
gap collapses from 4 eV at ao to essentially zero at 
the onset of the collapsed phase at w=0.55, mak- 
ing this an insulator-to-semimetal transition. The 
metallic phase then evolves continuously as the pres- 
sure is increased beyond P c . The 3d configuration 
can be characterized as t\ g = fc^ +*| s j.! resulting in 
a moment of lfis ■ The corresponding spin density is 
strongly anisotropic, although in a different manner 
than is the case for LDA+U. 



C. pseudo-SIC Method 

The spin-decomposed spectrum from this method 
is shown in Fig. |H1 symmetry-projected as done 
above for the LDA+U method. The PDOSs in the 
high-spin state are quite similar to those given by 
the LDA+U method. The transition could hardly 
be more different, however. The gap collapses in an 
insulator-to-good-metal character, the Fermi level 
lying within both majority and minority bands. The 
majority bands are the e' g pair (|m|=l) and are only 
slightly occupied. In the minority bands both e g and 
e' g are roughly quarter filled. The reason this solu- 
tion is (at least locally) stable seems clear: Ep falls 
in a deep valley in the minority DOS. 

In the collapsed, low-spin state, the a g orbital of 
both spins is occupied, and the majority e g pair is 
also fully occupied. This results in a configuration 
that can be characterized roughly as a g f e 2 g Vt a l I 
e"' 5 I (4) ' 5 I, giving spin § - § - \ = f . Thus the 
fact that the same moment is found in the low-spin 
state as was found with the LDA+U and pseudo-SIC 
methods seems accidental, because in those methods 
the energy gap required integer moment whereas the 
pseudo-SIC solution is firmly metallic. It is in fact 
close to half metallic, which accounts for the near- 
integer moment. In pseudo-SIC, the a g orbital is un- 
polarized (spin-paired), the e g pair {xz,yz} (|m| = l) 
is positively polarized, and the e' g pair {a; 2 — y 2 , xy} 
(\m\=2) is negatively polarized but to a smaller de- 
gree. 

D. SIC-LSD Method 

This method give much more tightly bound 3d 
state in the high-spin state than the other methods. 
An associated feature is that the majority-minority 
splitting, the "effective U" = U plus the exchange 
splitting 13 eV, more than twice as la large as used 
in the LDA+U method (both in this paper and else- 
where). Note that in the SIC-LSD method II s IC 
is a true Slater self-Coulomb integral, whereas in 
the LDA+U method the value of U represents the 
(somewhat screened) Coulomb interaction between 
a 3d electron and an additional 3d electron, so agree- 
ment between the two is not expected. Nevertheless 
the difference is striking. All five majority 3d states 
are localized, leading to the self-interaction poten- 
tial that binds them. The majority 3d states lie 6 
eV below the center of the 2p bands and hybridize 
very weakly, which accounts for the very narrow, al- 
most core-like 3d bands. 

In the collapsed phase, there are no localized 
states and the usual LSDA results re-emerges. All 3d 
states make some contribution to the moment, but 
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FIG. 8: (Color online) Spin- and symmetry-projected 
DOS of the high spin (upper panel) and low spin (lower 
panel) states resulting from the pseudo-SIC method. In 
the low-spin state, the a g orbital is unpolarized due to 
occupation by both spin directions, and shows little ex- 
change splitting. One of the e 9 -symmetry pairs, here 
called e' g , is unpolarized due to being unoccupied in both 
spin directions. It is the one labeled e g that is spin split 
across the Fermi level and is responsible for the moment. 



the strongest contribution arises from the a g orbital. 
One might argue in the LDA+U result that the net 
moment arises from the a g orbital (with the mo- 
ments arising from the e g and e' g orbital canceling). 
However, the electronic structure in the collapsed 
phase of this SIC-LSD method is very different from 
that of the LDA+U method. 



VII. DISCUSSION AND SUMMARY 

The results of four correlated band methods - 
LDA+U, SIC-LSD, pseudo-SIC, and HSE - have 
been compared for the equation of state (for both 
normal and collapsed phases), electronic structure 
(including 3d configuration), and the Mn moment 
under pressure. In order to make the comparison 
as straightforward as possible, the crystal structure 
was kept cubic (rocksalt). To compare seriously with 
experiment, one must account for the coupling of 
the AFM order to the structure, because this re- 
sults in a substantial rhombohedral distortion of the 
lattice. Then structural transitions (particularly 
to the B8 phase) must also be considered. 

The large volume, high-spin phases are qualita- 
tively the same for the various functionals: AFM 
with a fully-polarized 3d 5 configuration. Due to 
the large charge transfer energy, the configuration 
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FIG. 9: (Color online) Spin- and symmetry-projected 
DOS of the high spin (upper panel) and low spin (lower 
panel) states given by the SIC-LSD method. Orbital 
characters are expressed in the rhombohedral coordinate 
system. In the high-spin state, the majority 3d states are 
are centered 6 eV below the center of the occupied 2p 
bands, resulting in little hybridization and very narrow 
bands. The exchange splitting of the 3d states is about 
13 eV, providing an "effective U" from this method. 
The collapsed moment phase is representative of a band 
(LSD) ferromagnet. 

remains d° at all volumes studied here. The pre- 
dicted equations of state (which give the equilibrium 
volume, bulk modulus, and its pressure derivative) 
show rather strong variation, suggesting that the ex- 
tension of the Mn 3d and O 2p functions, or their 
hybridization, differ substantially even in this large 
volume phase. Of course, since the functionals are 
different, any given density would lead to different 
energies. 

Under pressure, the gap initially increases (all 
methods give this behavior) , and the system suffers a 
first-order transition (isostructural, by constraint) to 
a collapsed phase where hybridization must be cor- 
respondingly stronger. Uniformly among the meth- 
ods, the moment collapse reflects an S=| to S=i 
transition, rigorously so for the LDA+U method for 
which the collapsed phase retains a gap, and less rig- 
orously so for the the other methods where the col- 
lapsed phase is metallic. It is remarkable that none 
of the methods gives a collapse to a nonmagnetic 
state, which probably would be the most common 
expectation. 

This S=| — > o moment collapse is related in some 
cases to the local symmetry of the Mn 3d orbitals in 
the AFM phase (being most obvious for the LDA+U 
results of Fig. [fJJ). The symmetry is a g + e g + e' g , i.e. 
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a singlet and two doublets per spin direction. With- 
out further symmetry-breaking (orbital ordering) an 
S=| — > | transition requires a single spin flip, which 
could only be the a g spin. However, the a g state is 
more tightly bound than as least one of the two dou- 
blets both for LDA+U (Fig. © and for pseudo-SIC 
(Fig. |SJ. In LDA+U the e g doublet flips its spin, 
while in pseudo-SIC the a g singlet flips its spin leav- 
ing the minority e g and e' g doublets partially occu- 
pied and therefore metallic. This symmetry-related 
behavior depends of course on the magnetic order- 
ing that gives rise to the (electronic) rhombohedral 
symmetry. Above the Neel temperature, the mo- 
ment collapse at the Mott transition may proceed 
differently because the Mn moment would lie at a 
site of 'cubic' symmetry (a dynamic treatment could 
include the effect of short range spin correlations). 

Our study provides some of the first detailed in- 
formation on how magnetic moments in a real ma- 
terial may begin to disintegrate without vanishing 
identically, at or near a Mott transition, when cor- 
relation is taken into account. It is accepted that dy- 
namic processes will be required for a truly realistic 
picture of the Mott (insulator-to-metal) transition. 
However, a moment collapse between two insulating 
phases (as described here by two of the methods) 
may be described reasonably by a correlated band 
(static) approach. 

Beyond this similar amount of moment collapse, 
the four functionals give substantially different col- 
lapsed phases: differences in the Mn 3d magnetic 
configuration (although all remain d 5 ) and differ- 
ences in conducting versus insulating behavior. It 
is not surprising therefore that the collapsed-phase 
equations of state differ considerably between the 
methods. 

The differences in predictions can be traced, in 
principle, to the different ways in which exchange 
and correlation are corrected with respect to LDA. 
One clear shortcoming of LDA is in the local ap- 
proximation to the exchange energy. The HSE 
method deals with this problem directly, by using 
25% Hartree-Fock exchange. The self-interaction of 
the SIC-LSD method is largely a self-exchange en- 
ergy correction, subtracting out the spurious self- 



Coulomb energy that occurs in the Hartree func- 
tional if an orbital chooses to localize. A self- 
(local density) correlation correction is also included 
in SIC-LSD. The pseudo-SIC method includes the 
same correction if applied to an atom, but in a crys- 
tal the pseudo-SIC energy correction and change 
in potential takes a substantially different form, as 
the difference in predictions reflects. The LDA+U 
method is rather different in this respect: it specifi- 
cally does not subtract out any self-interaction (al- 
though it is sometimes discussed in this way). In 
the form Eq. of LDA+U correction, the sec- 
ond term is simply an LDA-like average of the first 
term. The on-site Coulomb repulsion is treated 
Hartree-Fock like, leading to an orbital-dependent, 
occupation-dependent potential. Each method has 
its own strengths, and each is only an anticipated 
improvement on LDA toward a better, more general 
functional. It is expected that more details of the re- 
sults may be published separately by the respective 
practitioners. 
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